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We study theoretically the feasibility of amplification and generation of terahertz radiation in 
dc-ac-driven semiconductor superlattices in the absence of electric domains. We find that if in ad- 
dition to dc bias a strong THz pump field is applied, Bloch gain profile for a small THz signal can 
be achieved under conditions of positive static differential conductivity. Here the positive differ- 
ential conductivity arises, similarly to the case of large-signal amplification scheme [H. Kroemer, 
cond-mat/0009311)], due to modifications of dc current density caused by the application of high- 
frequency ac field [K. Unterrainer et al., Phys. Rev. Lett. 76, 2973 (1996)]. Whereas the sign of 
absorption at low and zero frequencies is sensitive to the ac fields, the gain profile in the vicinity of 
gain maximum is robust. We suggest to use this ac-induced effect in a starter for THz Bloch oscil- 
lator. Our analysis demonstrates that the application of a short THz pulse to a superlattice allows 
to suppress the undesirable formation of electric domains and reach a sustained large-amplitude 
operation of the dc-biased Bloch oscillator. 
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I. INTRODUCTION 

Terahertz radiation (0.3—10 THz) has many promising 
applications in different areas of science and technology 
such as space astronomy, wideband communications and 
biosecuritji One of the main challenges in the develop- 
ment of THz technology is the construction of coherent, 
monochromatic and miniature sources of THz radiation 
that can operate at room temperature. A great recent 
achievement was the development of quantum cascade 
lasers that operate at THz frequencies employing quan- 
tum transitions between energy levels in multiple quan- 
tum well heterostructures 2 ". Continuous improvements in 
the design of THz quantum cascade lasers have allowed 
to increase their operational temperature above 100 K 3 . 
However, it is certainly very difficult to maintain popu- 
lation inversion in a THz quantum cascade laser at room 
temperature. 

Semiconductor superlattices (SLs) 4 , which are work- 
ing in the miniband transport regime, attract much at- 
tention as an artificial nonlinear medium demonstrating 
a resonant interaction with THz radiation (for review, 
see^&i). In the presence of dc bias, electrons from a 
single miniband of SL can perform transient THz Bloch 
oscillations&S, which decay at a picosecond time scale due 
to unavoidable scattering^. In the stationary transport 
regime, when the scattering of the miniband electrons is 
an important factor, the possibility of inversionless am- 
plification and generation of THz radiation in SLs has 
been well recognized after the seminal suggestion s 4 ' 11 . 
Bloch oscillations can potentially lead to the amplifica- 
tion of a weak ac field at frequencies smaller than the 
Bloch frequency, whereas absorption occurs for frequen- 
cies which are larger than the Bloch frequency^ 1 -. If elec- 
trons can perform relatively many cycles of Bloch oscil- 
lations between scattering events, the maximum of the 
small-signal gain is achieved in the vicinity of the Bloch 



frequency. This Bloch gain profile, shaped as a familiar 
dispersion curve, has recently attracted a lot of atten- 
tion because estimates predicted a significant THz gain 
near the Bloch frequency at room temperature. More- 
over, Bloch gain for THz signals of a large amplitude has 
also been estimate d 12 ! 13 . 

Within the semiclassical picture, the physical mech- 
anism of the high-frequency gain can be understood in 
terms of ballistic trajectories of electrons in the quasi- 
momentum space. It was found that the miniband 
electrons are gathered into certain favorable trajecto- 
ries forming electron bunches, which are eventually re- 
sponsible for the negative absorption 14 . Alternatively, 
the Bloch gain can be explained within theories which 
consider the scattering-assisted transitions between the 
Wannier-Stark states 1 ^. It is likely that the second order 
scattering-assisted mechanism of gain is not restricted to 
the case of single-band transport in SLs. The dispersive 
profile of optical gain has been recently observed in the 
quantum cascade lasers 1 ^. 

However, at least in the case of miniband transport in 
dc-biased SL, the small-signal gain profile extends to ar- 
bitrary low frequencies^ 1 -, which in essence means that 
the dc differential conductivity is also negative. SLs 
with static negative differential conductivity (NDC), sim- 
ilarly to the bulk semiconductors demonstrating Gunn 
effect, are unstable against the development of space- 
charge instabilit y 11 ! 18 . This electric instability eventu- 
ally results in the formation of moving domains of high 
electric field inside the SI>^. So far, the NDC-related 
instability has prevented an unambiguous observation of 
the net Bloch gain in long SLs. 

Importantly, it is still potentially possible to avoid the 
formation of the electric domains in SL, if the amplified 
THz signal has enough large amplitude. In this case, as 
has been demonstrated for a particular set of amplitude 
and frequency 2 ^, the dependence of the time-averaged 
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current on the applied bias can have a positive slope for 
a proper choice of the operation point. This effect in SLs 
is similar to the so-called Limited Space-charge Accumu- 
lation (LSA) mode of operation known in semiconductor 
devices with hot electron a 21 i 22 ' 23 . The finding 2 ^ can po- 
tentially be used for the construction of a THz amplifier 
with a large offset. However, to realize the THz oscillator 
it is still necessary to understand how to avoid the de- 
structive electric domains at small amplitude in order to 
reach the large-amplitude regime supporting stable oper- 
ation of the device. This is known as the device-turn-on 
problem for THz Bloch oscillator—. 

In recent years a number of interesting suggestions for 
the realization of small-signal Bloch gain in a SL were 
put forward. The first direction in this research area con- 
sists in the experimental investigations of the Bloch gain 
at short space or time scales, when the electric domains 
have not enough time to build up. Advanced nanostruc- 
ture design, an array of short SLs, has been introduced 
to measure a frequency dependent crossover from loss 
to gain 2 -. Modern ultrafast optical techniques were ap- 
plied to measure the Bloch gain in undoped superlattices 
during a short time window after a femtosecond optical 
excitation of carriers in SL 25-26 . Another interesting sug- 
gestion is to work with 2D structures, where the electric 
domains are effectively more suppressed in comparison 
with the case of 3D structure d 27 ' 28 . With this aim a lat- 
eral surface SL shunted by another SL^ have been intro- 
duced as an active medium for a future realization of the 
Bloch oscillator—. At least theoretically, it is also pos- 
sible to obtain a high-frequency gain in conditions of a 
positive dc differential conductivity (PDC) by engineer- 
ing the miniband dispersion relations in SLs^I or using a 
hot electron injection into a miniband^ 2 -. 

In this paper, we show theoretically how to reach gain 
for THz signals of both small and large amplitudes in a 
single miniband SL combining the action of constant and 
alternating electric fields. We focus on the requirements 
for gain under conditions of the electric stability. With 
this aim we first suppose that the electric field, acting on 
the miniband electrons, is Ed c + Eg coswi and re-examine 
the scheme of electrically stable large- amplitude (Eg) op- 
eration of Bloch oscillator 20 . We demonstrate that the 
gain in conditions of static PDC exists in rather wide 
ranges of w and Eg. In a real device, the ac probe field 
E$ cos tot should be the mode of a high-Q resonator. 

As a next step, we consider the action of the bichro- 
matic field Ed c + E ac cos fit + Eg cos u>t, where E ac cos fit 
is an additional external THz field (pump) and the am- 
plitude of the probe field (resonator mode), Eg, is now 
small. The pump Q and probe lj frequencies are assumed 
to be incommensurate. The basic idea behind this part 
of our paper is to use the external ac field to suppress the 
electric instability in SL for some range of bias Ed c . We 
show that in the case of dc-ac-driven SL the dispersive 
gain profile stays almost unaltered at high-frequencies 
lj ~ nfl (n = 1,2,...), although at small frequencies the 
absorption can change its sign and become positive (see 



Fig. [6]). Thus, THz gain without the electric instability 
can be achieved, and this bichromatic scheme can poten- 
tially be used to realize a THz generator. The expense 
we must pay for this new attractive possibility is the ne- 
cessity to use the pump THz field. 

We find that the physical mechanism of suppression of 
the electric instability in both these schemes is based on 
local modifications of the voltage-current (VI) character- 
istics of the SL induced by the THz fieldS2*2i. However, 
while in the large-amplitude mode of the Bloch oscilla- 
tor such changes in the dc current are self-induced by 
the strong resonator field Eg cos cot, in the bichromatic 
scheme, at least in the limit of weak probe field (Eg — > 0), 
they are mainly caused by the action of the external ac 
field. 

Our analysis also shows that it is potentially possi- 
ble to combine these two schemes of generation in a 
single device. Namely, we suggest to use the external 
field as a starter for THz Bloch oscillator: Temporal ac- 
tion of the ac pump can provide both the necessary sup- 
pression of electric domains and the gain for the small- 
amplitude oscillations of the resonator mode, before the 
field strength in this mode, Eg, can reach up to the am- 
plitude supporting electrically stable operation without 
external pump. Since the incommensurability of frequen- 
cies is supposed, the pump field of the starter should 
not necessarily be monochromatic. Intensive, broadband 
THz pulses (T-rays)2£ can be also used to ignite the sta- 
ble large-amplitude operation of the superlattice Bloch 
oscillator. Thus our theoretical research contributes to 
the solution of the device-turn-on problem of the canonic 
Bloch oscillator. 



II. NONLINEAR ELECTRON TRANSPORT 
AND THE CANONIC BLOCH OSCILLATOR 

Everywhere in this paper we work within the semiclas- 
sical approach based on the use of Boltzmann transport 
equation for the miniband with the tight-binding disper- 
sion relation^. We mainly use the standard approxima- 
tion of a single constant relaxation time r and employ an 
exact formal solution of the Boltzmann equation^; the 
effects of two distinct relaxation times are discussed in 
Appendix. We are interested in the electron dynamics 
under the action of a time-dependent field which consists 
of a constant and alternating parts E(t) = Ed c + E a i t (t). 
In general case, the ac part E a i t (t) can contain many 
incommensurate frequency components. We define the 
absorption of the probe ac field Eg cos Lot in the SL mini- 
band as 

^(w) = y(t)008Wt) tl (1) 

where j(t) is the current density induced in the SL by 
the total field E(t). In general case the averaging (. . .) f is 
performed over infinite time. Gain corresponds to A. < 0. 
In this section we will consider only the monochromatic 
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ac field 



E(t) = E dc + E s cosu>t. 



(2) 



For this field averaging in Eq. |1} should be performed 
only over the period T = 2tt juj. 

The absorption (gain) a in units cm -1 is related to A 
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where eo and c are the permittivity and the speed of light 
in vacuum, and n r is the refractive index of SL material. 
Whenever the gain overcomes the loss in the resonator, 
we have an oscillator. For the oscillator the generated 
power density V inside the sample is 

V = AE S . 

To estimate the gain and generated power in physical 
units we will use everywhere in this paper the following 
parameters of typical GaAs/AlAs SL: period d = 6 nm, 
miniband width A = 60 meV, density of electrons N = 
10 16 cm~ 3 , refractive index n r = \/l3 (GaAs). We con- 
sider operation at room temperature. The characteristic 
scattering time of miniband electrons is t ~ 200 fs. 

Before proceeding with the analysis of gain in ac-driven 
case, it is worth to remind the main nonlinear transport 
properties of a dc-biased SL. The dependence of the dc 
current density on the dc bias is given by the Esaki-Tsu 
formula^ 



jdc(eE dc d) = j p 
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2eE dc d/T 

i + (eE dc d/ry 
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Here V = 

E dc = E r . 



H/t and the peak current, corresponding to 
(E cr = h/{edr)), is 
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where vq = Ad/(2h) is the maximal electron velocity in 
the miniband and Ik(x) (k = 0, 1) are the modified Bessel 
functions. This temperature dependence of dc current is 
in a good agreement with the experimental. The depen- 
dence of j dc on E dc is shown in Fig. [1] Instead of the field 
strength variable E dc , the Esaki-Tsu characteristic can be 
represented in terms of Bloch frequency ujb = edE dc /H 
using the equality E dc /E cr = ujbt. For ujbt ~ 1 the 
Bloch frequency belongs to THz frequency ranged. 

The dc differential conductivity, a dc — dj dc /dE dc , de- 
fines the slope of VI characteristic at the working point. 
For the Esaki-Tsu dependence ([3]) it is 



a dc (eE dc d) 



2 jpeak 1 - (eE dc d) 2 /T 2 

[i + (eE dc d)yr 2 y 



(4) 



If E dc /E cr = u>bt > Ij °dc is negative and therefore 
Esaki-Tsu characteristic demonstrates NDC. 




FIG. 1: (Color online) The Esaki-Tsu characteristic. Figure 
demonstrates the quantum derivative for the frequency iv = 
3r _1 taken at the operation point E dc = 5E cr . Absorption at 
frequency u> is proportional to the difference quotient, Eq. (JS|, 
which is the slope of the inclined (red online) segment. Inset: 
Small-signal gain as a function of u for E dc = 5E cr . Bloch 



gain exists almost up to lo ~ ujb 
of gain at uj ~ ujb — t 1 = 4r _1 . 
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A. Small-signal gain and electric instability in dc 
biased superlattices 

We start with the consideration of small-signal gain 
in dc-biased superlattice @. In this case, the analytic 
calculation of absorption A(uj) gives the well-known for- 
mula of Ktitorov et alJ^. However, having in mind our 
future treatment of SL in bi- and polychromatic fields, 
it is instructive to introduce here the notion of quantum 
derivative and briefly analyze the high-frequency gain us- 
ing this tool. The absorption |T]) can be represented as 
the quantum derivative of the Esaki-Tsu current-field de- 
pendence © as&2&22 
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1 j dc (eE dc d + frhj) - j dc {eE dc d - huj) 

2 2hw 



eE s d. (5) 



Figure Q] demonstrates the meaning of quantum deriva- 
tive (finite difference) in the case of high-frequency probe 
field with foju = ST (ujt — 3) and the Wannier-Stark spac- 
ing Hub — 5r (E dc — 5E cr ). Geometrically, the quan- 
tum derivative represents the slope of segment with the 
length defined by the probe frequency to. The ends of 
the segment belong to the Esaki-Tsu curve and their lo- 
cations are determined by the choice of working point. In 
the quasistatic limit ujt C 1, we immediately get from 
Eq. (O the well-known result that the absorption is pro- 
portional to the dc differential conductivity at the opera- 
tion point (j4|) . On the other hand, as is obvious from the 
figure, in order to provide small-signal gain in the canonic 
Bloch oscillator the working point must be in the NDC 
portion of the Esaki-Tsu characteristic. 

The frequency dependence of absorption, calculated 
using Eqs. [5] and [3l is shown in the inset of Fig. [T] It 
illustrates the dispersive profile of Bloch gain at high 
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frequencies to ~ and also the existence of negative 
absorption for lot — ► 0. The frequency corresponding to 
the gain resonance can be most easily estimated for the 
choice of working point satisfying Ed c S> E cr . As is obvi- 
ous from the geometric meaning of Eq. ([5]), the maximum 
of gain is achieved when the left end of the red segment 
is located at the peak of VI characteristic. Then a simple 
geometrical analysis based on Fig. [T] immediately shows 
that the maximum of gain corresponds to the photon en- 
ergy hw w Hu>b — T {lot «wbt-1). This gain resonance 
indicates a dissipative quantum nature of the Bloch gain 
profile: Asymmetry in the elementary acts of emission 
and absorption is caused by the scattering. However, 
one should always remember that the feasibility of high- 
frequency generation in the canonic Bloch oscillator is 
complicated by the NDC-related electric instability. 

B. Photon-assisted peaks and large-signal gain 

In the case of large amplitude of the probe field, Eq.[2j 
the high-frequency gain is not anymore necessarily con- 
nected to the presence of static NDC, because the ac field 
can open up new transport channels leading to the for- 
mation of photon-assisted peaks in VI characteristi c 33 ' 34 . 
Here the dc current can be calculated with the help of 
photon replicas of the Esaki-Tsu current^ 

ti c (eE dc d) = Ji(0)3dc(eE dc d + Ihw) (6) 
i 

and the formula for absorption takes the form 

A = | ^ J, 09) [Ji+m + Ji-i(P)} j dc (eE dc d + Ihw), 
i 

(7) 

where J n (x) are the Bessel functions, summation is in 
infinite limits, /? = eEgd/fko and jd c is given by Eq. ([3]). 
We used the notation to distinguish the dc current 
modified by the action of ac field from the unmodified 
Esaki-Tsu current jdc- Using Eq. [5] it is easy to calculate 
the dc differential conductivity, tr^ c = dj% c j 'dEd c , 

°dc = J2 Jf{P)vdc{eE dc d + Ifrtu), (8) 
i 

where Ode is given by the formula ([¥]). 

Figure [2] shows the first photon-assisted peak in the 
dependence of j^ c on Ed c arising under the action of 
strong ac field. Importantly, the left wing of the peak 
(4 < E dc /E cr < 6) is characterized by PDC. PDC ex- 
ists if the Wannier-Stark spacing approximately equals 
to the photon energy \Hlob — fi^\ < T (other words, if 
\E dc / E cr - lot\ < 1). 

Now we want to demonstrate the feasibility of large- 
signal gain in conditions of static PDC. With this aim 
we choose the working point at the left wing of the 
first photon-assisted peak in Fig. [2] at Tilob = 5.5r 
{Edc = 5.5E cr ) and vary the amplitude of the ac field 
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FIG. 2: VI characteristic of a superlattice modified by the 
strong monochromatic field with Eg — 5E cr and lot — 5. 
Note the appearance of the photon-assisted peak centered at 
Edc/ E cr = 6. At this peak the differential conductivity is 
positive for 4 = lot — 1 < Edc/Ecr < lot +1 = 6. 

Eg. Absorption and, in our notations, power density will 
be negative if the SL can generate a high-frequency ra- 
diation. The dependence of the power density on Eg is 
depicted in Fig. [3] In the same figure we also indicated 
the range of Eg for which NDC is realized at the working 
point. We see that the generation definitely can be com- 
plicated by the existence of electric instability at small 
probe amplitudes: There always exists a threshold value, 
which the probe field amplitude must reach before SL can 
switch to PDC. However, if the small-signal space-charge 
instability can be suppressed by some way, for example 
by a very special SL design 24 , the generated power den- 
sity at frequency ui/2tt = 5t~ 1 /2it w 4 THz can reach 
~ 40 MW/cm 3 (Fig. [3]). For a typical semiconductor SL 
it corresponds to the power ~ 100 /xW. 

The whole ranges of bias Ed c and probe amplitudes 
Eg, supporting simultaneously PDC and gain at the 
first photon-assisted peak, are presented in Fig. [4] for 
one particular value of the photon energy fko = 5T 
{lot — 5). This figure shows that gain without NDC exist 
if Eg > 2E cr and hio < frio B < hio + T (5 < E dc /E cr < 6). 

Similarly, we also calculated the slope of jd c {Ed c ) de- 
pendence and the sign of absorption for many other probe 
frequencies. We summarize the requirements to obtain 
gain without NDC at the first photon-assisted peak as 

fuo > 2r, (9a) 
Eg > 2E cr , (9b) 
tiw < fuo B < hoo + T. (9c) 

The requirements (|9a|) and (|9bj) - which state that the 
photon energy must be larger than the scattering in- 
duced broadening and that ac field strength must be no 
less than two critical fields - guarantee the formation 
of a sharply defined photon-assisted peak. Note that in 
the limit opposite to Eq. (|9"aj) , the photon-assisted peaks 
never arise (see 4 ^ and references cited therein). Finally, 
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FIG. 3: (Color online) Generated power density as a function 
of the probe field amplitude Eg for fixed Edc = 5.5E cr and 
ujt — 5. The dark (red online) segment indicates the range 
of negative differential conductivity. If one assumes that the 
electric instability due to NDC (red online) at small probe 
amplitudes does not prevent the field in the resonator from 
growing, then the generation with positive differential con- 
ductivity can be achieved for Es E [2.5E C , W.8E C ]. 
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FIG. 4: (Color online) Regions of negative differential con- 
ductivity (red online) and gain (marked area bounded by 
the blue line) for ujt — 5. The region of gain overlaps 
with the region of positive differential conductivity (blank) 
for 5 = ujt < Edc I 'Ear < ujt + 1 = 6 and E 6 > 2E cr - The 
figure indicates the possibility to reach a large-signal amplifi- 
cation in dc-biased SL without NDC. 



the condition ([9c]) ensures that simultaneously (i) the ab- 
sorption is negative tuo < Hujb and (ii) the dc bias is 
chosen at the left wing of the photon-assisted peak cor- 
responding to PDC \tuvB— hu}\ < T. For completeness, we 
should also mention that PDC and gain can be simulta- 
neously achieved also by operating close to other photon- 
assisted peaks nfojj < Hojb < nfiw + T (n = 2, 3, . . .). In 
this case gain and modifications in the dc current den- 
sity occur due to multiphoton processes, and thus larger 
amplitudes of the probe field are required. 

We should note that a dynamical response of electrons 
to THz radiation can depend on details of the scatter- 



ing processes leading to the relaxation towards a thermal 
equilibrium. In the case of SLs with a single miniband, 
the simple model with two different scattering constants 
for the electron velocity j v and electron energy 7 e often 
becomes more adequate than the approach based on a 
single relaxation time^ 1 -. Our numerical calculations in- 
dicate that the inequalities © well describe the situation 
for ^e/lv > 0.5. However, a large-signal gain without 
NDC still can exist even if these scattering constants are 
very different (see Appendix [A")l . 



III. SUPERLATTICE IN BICHROMATIC FIELD 

We turn to the consideration of electron transport in 
a SL under the action of a bichromatic ac field 



E = Edc + E ac cos Sit + Es cos u>t. 



(10) 



We suppose that the ac pump E ac cosSlt is strong and 
the amplified probe field Es cos cut in general can have 
an arbitrary amplitude. In the case of generation, the 
probe field is a mode of a resonator. We assume that 
the frequencies SI and u> are incommensurate but both 
belong to THz frequency domain (Sir > 1, lot > 1). 

In our earlier discussion, the probe field of a large am- 
plitude induced some local structures into the VI char- 
acteristic of the SL, which eventually led to the large- 
signal gain without NDC. The main idea behind the 
bichromatic scheme of generation and amplification is 
to use an external ac field for the modification of VI 
characteristic* 3 ^!, thus making it possible to amplify 
even a small probe signal without NDC. We will also 
show that fields of large amplitude can be amplified in 
the same conditions as well. 



A. Main equations 

To begin with we need to calculate the dc current den- 
sity and absorption for the bichromatic field (jTDJ) of ar- 
bitrary amplitudes. By using the exact formal solution 
of Boltzmann transport equation 36 and following earlier 
contribution to the formalism 1 ^, we ge t 42 ' 43 for the dc 
current 



Jdc 



= J k(P)jdc(eE dc d + kfuv), (3 = eE s d/huj (11) 



and the absorption of the probe field 



k 



J fc+ i(/?) + Jk-i(P) 



j2 c (eE dc d + kfuv). 

(12) 

Here the dc current density modified by the pump field 
alone is 



jg(eE dc d) = J?(a)j dc {eE dc d + IKl), 



(13) 
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FIG. 5: (Color online) VI characteristic of SL under the action 
of a strong pump field (Sir = 5, E ac = 6E cr )- Additionally to 
the modified Esaki-Tsu peak, two photon-assisted peaks are 
visible. The operation point is chosen at the part of the sec- 
ond peak with a positive slope (Edc = 9.8E cr ). The quantum 
derivative (see Eq. ^15} ) at this operation point is negative 
for a wide range of ui values demonstrating the possibility of 
small-signal gain without NDC. Inset: Small-signal absorp- 
tion as a function of the probe frequency ui at Edc = 9.8E cr . 
There are two distinct gain resonances at frequencies uit ~ 3.8 
and lot ~ 8.8. 



where a — eE ac d/Ml. Next, the dc differential conduc- 
tivity ctJ" can be easily calculated from Eq. (fTTj) 

^dc = Ji(<x)J%(P)°dc(eE dc d + IMl + khu). (14) 

l,k 

NDC and PDC at the operation point Ed c correspond to 



a d c < and cr/" > 0, respectively. 

In the limit of small probe field (f) <C 1), we have from 
CD and ((111) that j% 



Jdc and 
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1 jSd eE dcd +_M - jg(eE dc d - hu) 

2 2hw 



eE s d, (15) 



where j£, is given by Eq. (flU]) . Therefore for a weak 
probe both dc current and absorption are determined 
by VI characteristic modified by the pump field. Since 
the formula for small-signal gain in the bichromatic case 
(fT5")) resembles the corresponding formula © after the 
substitution f d \ — ► jd c , we can directly extend the simple 
geometric analysis of sec. Ill Al to the case of more com- 
plicated modifications of VI characteristic caused by the 
action of strong ac pump. 



B. Small-signal analysis: Basic ideas 

We want to show how to achieve a high-frequency gain 
working at parts of modified VI characteristic with posi- 
tive slopes. With this aim we choose ac pump with large 
amplitude (E ac /E cr ^> 1) and high frequency (Sir » 1). 
For such pump fields new local structures in VI char- 
acteristic become most pronounced and simple analytic 



formulas describing the modifications of dc current and 
gain can be derived. Figure [5] illustrates the dependence 
of the dc current density on the dc bias which has been 
calculated using Eq. (fT5|) for Ml = hT (Sir = 5) and 
E ac — 6E cr . Additionally to the modified Esaki-Tsu 
peak one can see two new photon-assisted peaks. The 
peaks of the current are centered at 



Edc/E 

c 



1 



(16) 



with n = 0,1,2. The left wings of the photon-assisted 
peaks (4 < Ed c /E cr < 6 and 9 < Ed c /E cr < 11) have 
positive slope, i.e. PDC. We found that PDC arises if the 
Wannier-Stark spacing approximately equals to the en- 
ergy of one photon \Hujb — < T (Sir — 1 < Ed c /E cr < 
Sir + 1) or to the energy of two photons \ fouJB — 2Ml\ < T 
(2S1t — 1 < Edc/E cr < 2S1t +1). In a more general case 
the necessary condition for PDC becomes 



\Hlub -nW\ <T (n = 0,1,2...). 



(17) 



The choice of the operation point at these parts of VI 
should prevent the development of space-charge instabil- 
ity typical for systems with NDC. For example, in Fig. [5] 
we choose the working point at the left part of second 
photon-assisted peak. 

Following Eq. (fTSf the calculation of high-frequency 
absorption of a weak field requires the finding of quantum 
derivative of j dc at the working point Ed c - The geomet- 
ric meaning of the quantum derivative, which is evident 
from Fig. [51 is similar to the one described in sec. Ill Al 
Importantly, the slope of (red online) segment, which is 
determining the quantum derivative, can be negative for 
many to providing gain for these frequencies (Fig. [5j . 

It is easy to see that gain has a maximum if the left 
end of the segment is located at one of the dc current 
peaks. Taking into account that for high frequencies the 
segment has a small slope and the peaks are centered at 
Edc given by Eq. (|TB)) . we get that the gain resonances 
exist for the values of bias and probe frequency satisfying 
Edc/ E cr — ~ £;S!t + 1 (k = 0, 1, 2). In terms of photon 
energies and Wannier-Stark spacing, this condition for 
the gain resonances takes the form 



fkrf ~ hujB — kMl — r. 



(18) 



The inset of Fig. O shows the absorption as a function 
of uj for huj B = 9.8T, Ml = 5T and E ac = 6E cr . The 
gain profile demonstrates two distinct gain resonances at 
hui w 3.8r and hw as 8.8T, which are well described by 
Eq. CEHD with k = 1 and k = 0. 

The dispersive gain profile and the locations of gain 
maxima indicate that from a pure quantum perspec- 
tive, the gain in the dc-ac-driven SL originates from 
scattering-assisted transitions between quantum mechan- 
ical states. In contrast to the Bloch gain profile in a pure 
dc-biased SL, the photon sidebands also play their role 
here. 

An interesting question is how the dispersive shape of 
gain profile in dc-ac-driven SL is sensitive to the am- 
plitude and frequency of the pump field. Within the 
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FIG. 6: (Color online) Small-signal absorption as a function 
of the probe frequency u) for Qt = 5 and different values 
of the pump amplitude E ac /E C r =: 0,1.2,2.1,3. The gain 
profile near the maximum of gain changes only slightly with 
an increase of the pump amplitude. However, the absorption 
at low frequencies changes from negative to positive indicating 
the transition from NDC (online red curves) to PDC (online 
black curves) at the operation point Edc = A.c\E cr . 

geometrical representation of the quantum derivative a 
crossover from gain to absorption corresponds to the 
movement of the end of the virtual segment from peak 
to dip in VI characteristic. Therefore one can expect 
that the dispersive shape of the absorption profile near 
this crossover is universal as long as well pronounced lo- 
cal structures in VI do exist. Our numerical calculations 
indeed demonstrate that in the vicinity of the crossover 
point, which also includes the gain resonance, the gain 
profile remains almost unaltered. 

Figure [5] shows the small-signal absorption as a func- 
tion of uj for different amplitudes of the pump, E ac 6 
[0, 3-E cr ], with the high frequency hil = 5T. We choose 
the operation point E& c — A.8E cr . If the ac pump 
is strong enough, this operation point belongs to the 
first photon-assisted peak with PDC (see Fig. [5l where 
E ac = &E cr ). For all amplitudes the dispersive Bloch 
gain profile at high-frequencies is maintained together 
with with the gain resonance at the photon energy hw rs 
Hujb — r = 3.8r. However, with increasing pump ampli- 
tude the low-frequency absorption changes from negative 
to positive; the switch from NDC to PDC takes place at 
E ac w 2E cr . Although the electric stability is achieved by 
a drastic change in the low-frequency absorption, there 
is practically no expense in the high-frequency gain. 

To have an analytic formula describing the require- 
ments for the existence of Bloch gain at PDC in dc-ac- 
driven SL, we directly combine the condition for PDC 
(fT7| with the formula for gain resonances (fPB)) . and get 
an inequality (n — k)Ml — 2T < fku — Hujb — kMl — V < 
(n — k)Ml. It can be rewritten as 

\huj -mhn + T\ <T, (19) 

where the index m (m = 1,2,3...) marks the m-th range 



of gain in the absorption profile. We observed numeri- 
cally that this naive estimate works surprisingly well for 
high-frequency pump fields with amplitudes E ac > 2E cr . 
Here the underlying physics is the robustness of gain pro- 
file at high frequencies with its simultaneous sensitivity 
to variations of the pump field at low frequencies. 



C. General conditions for gain: signals of arbitrary 
amplitude and positive dc conductivity 

We want to show how to extend our analysis to the 
case of probe field which is not weak anymore. Our com- 
putations demonstrate that for parameters satisfying the 
small-signal gain resonance, the large-signal gain as a rule 
also exists. Moreover, we found that the gain profile near 
gain resonances stays almost unaltered as the probe am- 
plitude changes. Furthermore, the important condition 
(fT9|) works also at large amplitudes of the probe field. 
However, simultaneously the analog of Eq. (fTT)) for the 
probe field must be also introduced: Wannier-Stark spac- 
ing is approximately equal to an integer multiple of the 
photon energies at the probe frequency, \fno b — lfno\ <Y. 

Now we can summarize our findings and list all nec- 
essary requirements to be satisfied in order to use the 
bichromatic scheme for a generation of high-frequency 
field under conditions of electric stability. 





m 


> 


2T, 




(20a) 




E ac 


> 


2E C r , 




(20b) 




-r < 


Hujb 


< nMl -\ 


-r, 


(20c) 




hu> 


> 


2r, 




(20d) 


IhjJ 


-r < 


fojJ B 


< Ihjj + 


r, 


(20e) 




2r < 




< mhfl, 




(20f) 



where in general integers n, I, and m are different. Some 
of these inequalities, like for example (|20b|) and (|20cp . 
are of course same as derived earlier and here are listed 
for completeness. We should note these conditions are 
asymptotic and work better the larger is the amplitude 
(Eq. I20b[) and the higher are the frequencies (Eqs. I20a[ 

EM]). 

The requirements (|20all20c[) are necessary for the am- 
plification of a weak probe field at an incommensurate 
frequency ui in conditions of PDC. They ensure that the 
operation point is located at the part of well-defined n- 
th photon-assisted peak having a positive slope. In order 
the probe field still can be amplified even if its amplitude 
is not small anymore, the requirements (|20dll20fj) should 
be taken into an account. Below we will demonstrate how 
the requirements (|20|) work in the cases of THz genera- 
tion at the first and second photon- assisted peaks. The 
effect of different relaxation constants, j v ^ 7 e , is ana- 
lyzed in Appendix [XJ 
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FIG. 7: (Color online) Regions of NDC (red online) and gain 
(marked, bounded by blue curve) for E ac = 2E cr , Qt — 5.5, 
and lot = 5. The figure illustrates the feasibility of gain 
in conditions of suppressed electric instability for the wide 
ranges of signal amplitudes Es and bias values Edc under the 
action of moderate ac pump. 



1. Gain and generated power density at the first 
photon- assisted peak 
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Figure [7] shows the regions of gain and NDC for dif- 
ferent values of dc bias and probe amplitude Eg, when 
the pump amplitude is fixed (cf. Fig. 0J. It demon- 
strates the possibility to achieve simultaneously PDC 
and gain at the first photon-assisted peak [n = 1 in 
Eqs. ([20]) ] for the probe frequency in the proximity of 
the pump frequency [m = 1 in Eqs. (|2"0)) ]. In this fig- 
ure we choose a relatively low amplitude of the pump 
E ac = 2E cr , which is close to the threshold amplitude 
[see Eq. (|20b[) ]. As a consequence, PDC for E$ ^ ex- 
ists for 5.3E cr < Edc < 5.8E cr , what is a bit narrower 
than the range 4.5E cr < Edc < 6.5-E cr following from 
Eq. (|20c[) . On the other hand, figure [5£i shows the re- 
gions of gain and NDC in the plane E ac -Es for a fixed dc 
bias. It demonstrates that gain in conditions of PDC ex- 
ists for the pump amplitudes E ac € [2E cr , 7.6E cr ], when 
requirements of Eqs. (|2"D|) are satisfied. 

Figure [9] shows that the generated power density at 
the first photon-assisted peak in the bichromatic scheme 
slightly decreases with an increase in the pump ampli- 
tude but still remains comparable with the power den- 
sity which can potentially be generated in the canonic 
Bloch oscillator (cf. Fig. [3]). Moreover, a comparison of 
Fig. [5] and Fig. [IJ inset) allows to make similar statement 
concerning the magnitudes of gain in these two schemes. 



2. Generation with frequency shifting 

By applying a very strong pump field it is possible, 
due to multiphoton processes, to exploit other photon- 
assisted peaks and reach the generation at frequencies 
significantly different from the pump frequency. In the 
case of second photon-assisted peak (n = 2 in Eq. (|20cp ). 



FIG. 8: (Color online) Regions of NDC (red online) and gain 
(marked, bounded by blue curve) in the plane of pump-probe 
amplitudes for (a) Qt — 5, lobt = 4.8 and lot = 4 and (b) 
Qt = 5, ujbt = 9.8 and ujt — 9. Gain can exists without 
NDC. The arrows schematically illustrate the use of ac pump 
in the THz starter for Bloch oscillator (see sec. HVI for details). 




2 4 6 8 10 12 

8 cr 

FIG. 9: (Color online) Generated power density as a function 
of the probe amplitude for several increasing values of the 
pump. The black curves correspond to the generation at the 
first photon-assisted peak (Qt = 5, Edc = 4:.8E cr , lot — 4) 
induced by the ac pump with E ac /E cr =: 2.1, 3.5. The brown 
curves demonstrate that the generated power at the second 
photon-assisted peak (Qt = 5, Edc = 9.8i? C r, lot — 9) more 
strongly depends on the pump amplitude (E a c/E cr =: 6,7). 
As can be seen from Fig. [8] the generation in both cases can 
be realized without the electric instability. 
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the frequency of the probe field can be easily shifted 
from the pump frequency so that the photon energy will 
be close to the gain resonance hui w 2M7 — T, corre- 
sponding to m = 2 in (|20fj) . On the other hand, the 
requirement (|20e|) can be satisfied for I = 1 by tuning dc 
bias. Therefore, for enough high frequency of the pump 
all requirements (f20|) are satisfied, and amplification of 
the probe field with both small and large amplitudes can 
be achieved in conditions of PDC. Figures [3 [SJd, and [9] 
demonstrate this kind of possibilities for the pump fre- 
quency Or = 5 and dc bias Edc = 9-8E cr . 



IV. THZ STARTER FOR BLOCH OSCILLATOR 

It is easy to notice that the large-signal gain in con- 
ditions of PDC is realized at the same parameters of 
probe fields in both the monochromatic and bichromatic 
schemes of amplification. Good illustration of this comes 
from the comparison of figures [4] and [7] They show the 
regions of gain and PDC in these two different schemes 
for the same photon energy of the probe field fko = 5T. 
In both cases the regions of gain and PDC overlap ap- 
proximately for dc bias values 5.1E cr < Ed c < 5.9E cr 
and probe amplitudes 1.9.E cr < Eg < 11.5E cr . How- 
ever, in the bichromatic case there also exists an ad- 
ditional overlap region at small probe amplitudes for 
5.3E cr < Edc < 5.8E cr due to the formation of photon- 
assisted peak induced by the action of ac pump. 

In more general case, we can compare the requirements 
© and (j2"0|) . Assuming that the probe frequency and dc 
bias satisfy conditions (|9a[) and (|9c|) , we immediately see 
that the conditions (120dj) and (I20el with 1 = 1) are also 
satisfied. Moreover, the rest of requirements (j20|) can be 
satisfied with a proper choice of the pump field. 

Therefore it is potentially possible to combine both 
these schemes in a single device. The process of domain- 
free generation in a resonator with very high-Q will in- 
clude two stages: In the first stage, the ac pump is used 
to excite the field mode with an amplitude Eg > 2E cr 
(requirement (I9aj) ). In the second stage, the pump field 
is gradually switched off and sustained generation is per- 
formed at a large amplitude of E$. We underline that 
here the whole switching off process takes place in condi- 
tions of strong gain and PDC. Figure [5] illustrates these 
two stages of generation utilizing the first (a) and sec- 
ond (b) photon-assisted peaks. Arrows show a safe pass 
around the dangerous NDC region in the plane E ac -Eg, 
which corresponds to the ignition of canonic Bloch oscil- 
lator. 

Since the pump ac field is needed only to the initial 
stage in order to suppress the formation of electric do- 
mains, and it is then switched off, it seems reasonable to 
use THz pulse instead of CW field in this THz starter for 
Bloch oscillator. The required high-power pulses of 100 
ns - 1 fis duration can be generated, for instance, with 
the help of free-electron laser or C02-laser. 

A quite different situation arises when few-cycle THz 



pulses are in use. The field cannot be considered as a 
monochromatic anymore. Therefore now we will ana- 
lyze the feasibility of gain in the case of polychromatic 
ac pump. For simplicity, we consider the case of trichro- 
matic field 

E ac cos fit E ac cos^M 

pump — 2 4 

| E ac cosn 2 t 

+ 4 ' ^ Zi > 

The amplified probe field continues to be monochromatic 
Eg cos tot. We suppose that the pump Q, fii, £l 2 and 
probe, ui, frequencies satisfy the generalized incommen- 
surability condition. Namely, we assume that there does 
not exist any nontrivial integers m, (i — 1,2,3,4) such 
that mif2i + m 2 Q 2 + + m^io = 0. The equation 

(|21[) can be considered as a crude approximation for the 
temporal shape of a real broadband THz pulse, where 
the Fourier transform of the pulse has a maximum at f2 
and |Oa — is related to the bandwidth of the pulse. 

The analysis of previous sections can be directly gener- 
alized to the case of polychromatic pump field (f2"Tj) . and 
we get for dc current density 

JT 1 ^ = E J kW)jdc in2 i eE dcd + khui) (22) 

k 

and for absorption of the probe field 

A = lJ2 J ^){ J k+i(^)+Jk-im 

k 

j^ n2 (eE dc d + khw). 

(23) 

Here j^ c 0l ° 2 is dc current density in SL modified by the 
pump field alone 

j2o lQ2 = E J?^/2)J 2 m (a/4)J^(a/4) 

jdc(eE dc d + IMl + mteli + nhfl 2 ). (24) 

Equations (|2"2")l and (j2"3")) are valid for an arbitrary ampli- 
tude of the probe field. Remarkably, in the limit of weak 
probe field (Eg — > 0) , the absorption A becomes propor- 
tional to the quantum derivative of the pump-modified 
dc current (pH)) . 

In the limit Qj — ► f2 (i = 1, 2) the trichromatic pump 
field (|2"Tj) becomes momochromatic and therefore one 
would expect a restoration of all results of section IIIII 
However that is not the case, and the current density 
and absorption, given by the equations (|22|) and (|23[) . 
are not reduced to the corresponding equations of sec- 
tion [TTTJ The disagreement arises because the condition 
of incommensurability is not satisfied in this limit. 

We computed the regions of NDC and gain for the 
pump (|2ip with two different side frequencies f2i and 
when the cental frequency f2 determines the first photon- 
assisted peak in VI characteristic. Figure [TU] shows the 
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FIG. 10: (Color online) Regions of NDC (red online) and 
gain (marked) for lubt — 4.8, lot = 4 and the trichromatic 
pump field (|21[) with Sir = 5, Q\t = 5.5 and Q.2T = 4.75. 
Comparison of this figure with Fig. [S^, demonstrates that the 
realization of gain in conditions of PDC under the polychro- 
matic pump requires larger values of E ac than in the case of 
momochromatic pump. 



regions of gain and PDC in the E ac -Es plane. The 
Wannier-Stark spacing and photon energies (Hub = 4.817, 
hfl = 5T and hio = AT) are the same as used in figure [8^,. 
By comparing figures [5] and [TUl one observes that the 
gain and PDC can be achieved in both cases at similar 
amplitudes of the probe field, while a bit larger pump is 
required in the trichromatic case. 

These results support the possibility to replace a 
monochromatic pump with a broadband THz pulse in 
the starter for Bloch oscillator. Required intensive T- 
rays can be produced by a number of ways, including the 
use of interdigitated photoconducting device^ 4 -, amplifier 
systems^ 5 - or from ultrafast ionizing air— (for recent re- 
view, see 3 - 5 -). 



V. CONCLUSION 

We theoretically analyzed the feasibility to reach THz 
gain in dc-biased semiconductor superlattice at room 
temperature using mono- and polychromatic alternating 
pump fields. We showed that different kinds of pump 
fields can suppress the formation of high-field electric do- 
mains inside the superlattice while still preserve a broad- 
band THz gain at frequencies incommensurate with the 
pump frequencies. 

Our approach is based on the utilization of well-known 
photon-assisted peak o 33 i 34 arising in the voltage-current 
characteristics of THz-driven superlattices. Choice of the 
operation point at the positive slope of such peak al- 
lows to suppress the electric instability. For the search of 
high-frequency gain in conditions of the positive differ- 
ential conductivity we employed simple but powerful ge- 
ometric interpretations of the intraminiband absorption 
formulas. Combining these analytic tools and numerical 



simulations, we demonstrated that the Bloch gain pro- 
file in dc-ac-driven superlattice is robust near the gain 
resonances. 

We suggested to use the robustness property of THz 
gain in an effective starter for the canonic (only dc- 
biased) Bloch oscillator operating in the electrically sta- 
ble large-amplitude mode of generation^ . Since THz 
gain in both the canonic Bloch oscillator and dc-ac- 
pumped Bloch oscillator occurs in the same ranges of dc 
bias, generation frequency and amplitude, a temporal ap- 
plication of the ac pump field can be used for domain-free 
generation up to the amplitude, which becomes sufficient 
for a stable operation of the canonic Bloch oscillator. The 
pump field, which ignites the Bloch oscillator, can also 
be polychromatic. We predicted that available intensive, 
broadband THz pulses (T-rays) can potentially be used 
in the fast ignition of Bloch oscillator. 

We conclude with two remarks. First, in this work we 
limit our attention to the case of incommensurate fre- 
quencies. If the pump and probe frequencies are com- 
mensurate, an additional coherent term in the formulas 
for high- frequency absorption will arise 47 . The sign of 
the term depends on the value of relative phase between 
the pump and amplified neld a 47 i 48 . However, at least 
in the case of generation with the monochromatic pump 
field, the coherent term always contributes to the high- 
frequency gain. A more detailed analysis of this situa- 
tion, including a consideration of microwave and sub-THz 
pump field o 49 ' 50 , will be presented elsewhere. 

Second, we have focused on the local modification 
of voltage-current characteristics of superlattices arising 
when the Bloch frequency approaches the ac field fre- 
quency or its harmonics. However, similar modifications 
of the electric characteristics have been observed in the 
semiconductor superlattices subject to the tilted high 
magnetic field, when the Bloch frequency is close to the 
cyclotron frequency or its harmonics 51 . It is very inter- 
esting problem to consider modifications of the voltage- 
current characteristics 5 ^ and analyze the feasibility of 
THz gain in a superlattice under the combined action 
of tilted magnetic and alternating electric fields. 
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FIG. 11: (Color online) Regions of NDC (red online) and gain 
(marked) in dc-biased superlattice (Eq. ((2j) for the probe field 
with the frequency lot = 5 and for several values of v: (a) 
v = 1, (b) v — 0.4 and (c) v — 0.1. The plots illustrate the 
feasibility of large-signal amplification under suppressed elec- 
tric instability for various ratios of the scattering constants. 
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FIG. 12: (Color online) Regions of NDC (red online) and gain 
(marked) in ac-pumped superlattice (Eq. (JTOj) ) for fir = 5, 
lot = 4, ujbt = 4.8 and for several values of v: (a) v = 1, (b) 
v = 0.4 and (c) v — 0.1. The plots illustrate the feasibility 
of generation under suppressed electric instability for various 
ratios of the scattering constants. 



APPENDIX A: EFFECT OF ELASTIC 
SCATTERING 

Here we consider the effects of different scattering con- 
stants in both the monochromatic ([2]) and bichromatic 
(fTU|) schemes of THz amplification and generation. The 
following SL balance equations can be derived from the 
Boltzmann transport equation for a single miniband of 

d , , d 2 e 2 E(t) . . 
J t Ji f )+ h 2 — g (*) = -7vj(*). 

±s(t)-E(t)j(t) = - le [s(t)-e eq \, (Al) 

where j{t) is the current density and e{t) is the total 
energy density of electrons within the miniband, 7 £ and 



7u = 7 £ + 7 e i are the phenomenological scattering con- 
stants for electron energy and miniband electron veloc- 
ity 7 e i is the scattering constant describing elastic scat- 
tering events. All electrons are at the bottom of mini- 
band for e = —TV A/2 and the upper edge of miniband 
is reached if e — +NA/2. The average electron en- 
ergy in thermal equilibrium e oq depends on the lattice 
temperature^. It is convenient to introduce the mean 
scattering time r = l/^/7u7 e and ratio of scattering con- 
stants v — 7 e /7„ < 1. According to many experiments 53 
v > 0.5 is a good assumption, although it is not valid for 
all superlattices^. 

Solving numerically the SL balance equations (|A1[) . we 
determine the time-dependence of a steady current for 
t ^> t and then calculate the dc differential conductivity 
and absorption of the probe field. In the case v = 1, 
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we find an excellent agreement with the earlier analytic 
results following from the formulas (J7J), (JSJ), and (fT2")l . 
(fFfjl . We observe that the photon-assisted peaks occur 
in VI characteristics for all values of v > 0.1; however, 
with a decrease in the ratio v these VI structures become 
less pronounced and magnitudes of high-frequency gain 
are decreasing. 

Importantly, the amplification still occurs in approxi- 
mately the same ranges of field amplitudes for different 
values of v. Figure QT] represents the ranges of dc bias 
Edc and ac pump amplitude E ac , supporting large-signal 
gain with PDC at the first photon-assisted peak induced 



by the monochromatic field ([2]). With decreasing v the 
range of amplitudes for negative absorption shrinks but 
the range of amplitudes supporting PDC expands (cf. 
subplots in Fig. [TTj) . The changes are obviously small as 
long as v > 0.4. Figure [12] shows the ranges of pump 
E ac and probe Eg field amplitudes supporting a genera- 
tion at the first photon-assisted peak in the bichromatic 
scheme. The suppression of electric domains and the ab- 
sorption in the bichromatic field (flC)]) are based on the 
same physical phenomena. Therefore the decrease in v 
affects the bichromatic scheme in a similar way as in the 
mochromatic case. 
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